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We consider a model for electron tunneling between a pair of ultrasmall metallic grains. Under 
appropriate circumstances, non-equilibrium final state effects can strongly enhance tunneling and 
produce electron-hole coherence between the grains. The model displays a quantum phase transition 
between a Coulomb blockaded state to a coherent state exhibiting subohmic tunneling conductance. 
The critical state of the junction exhibits a temperature independent resistance of order h/e 2 . Finally 
we discuss the possible relevance to granular materials and quantum dots. In particular, similarities 
between the quantum transition in our model and the metal-insulator transition in granular wires 
observed by Herzog et al. are described in detail. 
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Recent work by Herzog et al. jjj have found dramatic 
evidence for an unusual metal insulator(MI) transition in 
granular wires fabricated by in situ deposition through 
a metallic stencil onto a GaAs substrate. The transition 
involves an abrupt multi-order change in wire resistance 
as a function of the amount of deposited metal. The 
transition occurs in a variety of materials (including Sn, 
Pb, Au, Ag, and Pbo.85Bio.15). Although the transition 
may be observed in wires as wide as 7000A, the resistance 
gap decreases with increasing wire width and is absent 
in two-dimensional films. 

In this paper, we wish to consider the possibility of 
modeling the above behavior in terms of a pair of metal- 
lic grains in close proximity. Our model will include 
two ingredients. The first is electrostatic charging ef- 
fects. This should be important in the Herzog et al. 
experiments since the charging energy Eq is estimated 
|^| to be large (100K). The second ingredient is a non- 
equilibrium final state effect in which electrostatic fields 
between the two particles are suddenly switched on dur- 
ing the tunneling process. This effect is an exciton ef- 
fect in which the tunneled electron is attracted to the 
positively charged counter-electrode. We will show that 
the competition between the exciton effect and Coulomb 
blockade gives rise to a MI transition between a phase 
exhibiting sub-Ohmic I(V) characteristics to a phase ex- 
hibiting a Coulomb blockade. The critical state separat- 
ing these two phases exhibits a temperature independent 
conductance. Finally, we will discuss the the similarities 
between these results and those of Herzog et al. Q . 

We begin our discussion by considering a pair of identi- 
cal metallic grains on an insulating substrate (see Fig. [I]). 
The grains are in close proximity and form an ultrasmall 
tunnel junction with intergrain capacitance C < 10~ 15 F 
(For the present argument, wc will neglect the intragrain 
capacitance). The two grains may be part of a granular 
host. However, tunneling to other grains in the host ma- 



terial will be ignored. Now consider a tunneling process 
in which an electron tunnels from grain #1 to grain #2 
(see Fig. lc-ld). In the classical Coulomb blockade pic- 
ture a tunneled electron causes all energy levels in the 
grain #2 to up-shift by e 2 /2C and all energies in grain 
#1 to down shift by e 2 /2C. In the diagram, we have 
represented the electrostatic potentials and surface con- 
finement potential as a square well whose shape is unaf- 
fected by the tunneling process. Notice that the classical 
Coulomb blockade picture does not properly describe the 
non-equilibrium effect associated with suddenly switch- 
ing on the electrostatic attraction between the two grains 
(see Fig. le-lf). Such effects can give rise to shake-up 
(orthogonality catastrophe) effects which can seriously 
effect tunneling rates especially in small particles. 
To see how important these effects might be, consider 
the classic problem of the x-ray absorption edge (||| . In 
that problem, the absorption intensity, near threshold, 
could vanish at threshold due to the orthogonality catas- 
trophe H or could exhibit a power law divergence known 
as the exciton effect A priori one might expect that 
a similar effect could cause the differential conductance 
of a tunnel junction to vanish or diverge. 

To examine these effects in detail, we consider a 
model which describes tunneling between the 

two grains, 
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where Ht — Sfc P Pfcp e c fc c p + n - c -l an< ^ where <f> is de- 
fined in terms of the voltage difference Vl — Vr across 
the grains via the relation </> = c(Vl — Vr)/Ti. Finally, 
the Hl and Hr are given by 
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FIG. 1. (a) An illustration of two nearby metallic grains 
on an insulating substrate. The two grains form an ultrasmall 
tunnel junction. Although the grains may be part of a com- 
posite granular material, coupling to other grains in the host 
material will be ignored, (b) Energy levels and confinement 
potential associated with (a) . (c) Metallic grains after tunnel- 
ing event. An unrealistic charge distribution is obtained when 
V kk i is set to zero, (d) The energy levels and confinement po- 
tentials associated with (d). (e) Metallic grains after tunnel- 
ing event. Charges localized at the tunneling site give rise to 
long range electrostatic interactions which act as a suddenly 
switched-on potential V kk i = V£.t. (f) The non-equilibrium 
energy levels and confinement potential associated with (e). 

where V kk , represents the sudden change of surface and 
electrostatic potentials on grain a — L, R which occurs 
during the tunneling process. If V kk , — 0, we recover 
the standard model (see Figs, lc-ld and 2a) which 
describes the effects of particle hole excitations induced 
by tunneling processes in ultrasmall tunnel junctions. 
V kk , — V k \, ^ would be chosen to include shake-up and 
other final state effects in symmetric tunnel junctions. 

Now consider the zero temperature tunnel conductance 
to leading order in T kp — T. One may calculate the 
tunneling current using / = —2elm{X ret (—eV)} where 
X ret = \T\ 2 cItG r (t)G , 1(t) expiajr is the retarded 
response function associated with the non-equilibrium 
production of the R-electron and L-hole. Now the 
non-equilibrium electron and hole propagators associ- 



FIG. 2. Differential conductance for e < (dot-dashed), 
e = (solid) and e > (dashed), e > curve assumes exciton 
effect dominates the orthogonality catastrophe. 

ated with the suddenly switched on V kk , are given by 

G h jl e {t) oc N R t-( 1±SR W exp -iEft where N a is the 
DOS at the Fermi level in the a = L,R electrode and 
S a is a phase shift associated with the scattering of elec- 
trons off the potential V kk , . In the above result, + is used 
for hole and — for electrons in the right electrode. Simi- 
larly G h L ' e {t) oc N L t^ SL /^ 2 exp -iEft for the left elec- 
trode. Because of intergrain charging energy, the elec- 
trons tunnel to a non-equilibrium state characterized by 
Fermi levels shifted such that Ef — Ej = e 2 /2C. Com- 
bining these results one obtains ^ oc ^ylji^c ) 

where e = 2(8r/tt + Sl/k) ~ (8 R /ir) 2 — (8l/^) 2 and 
R^, 1 — Aire 2 NlN r \T\ 2 For small positive 5 a , there is 
a competition between excitonic effects associated with 
2(8r/tt + 8l/'^) and orthogonality effects associated with 
— (Sr/it) 2 — (Sl/tt) 2 . Depending on which terms dom- 
inate, one can obtain (Fig. 2b) a divergent differential 
conductance || at threshold Vr = e 2 /2C or (Fig. 2c) a 
vanishing conductance at threshold. 

In general the form of the final state interactions V kk , 
is not important to the following discussion and is diffi- 
cult to calculate. However there are several observations 
which can be made: First we observe that one can esti- 
mate p0| , pj| \5l,r\ ~ it/k where k is the number of trans- 
verse channels available for intergrain tunneling. Next, 
we observed that 6r will be positive if V kk , is a potential 
which tends to keep the tunneled electron in the electrode 
near the tunneling site i.e. if the electron is attracted 
to the positively charged electrode. Similarly Sl will be 
positive if the hole is attracted to the negatively charged 
electrode. Hence the electrostatic interaction between a 
pair of grains with a small number of accessible tunneling 
channels is expected to give a positive 5 a large enough 
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to make exciton effects observable. 

To some readers it may be surprising that the repulsive 
Coulomb interactions would enhance tunneling between 
the grains. The behavior is not unusual and can be found 
in several simple models. For instance, consider a pair of 
semi-infinite 1-D spinless chains described by the Hamil- 
tonian H = [t'c^cm + h.c] + J2i=ol tc li+i c Li + ^- c -] + 
C/(c+(0)c L (0) - l/2)(4(0)c fl (0) - 1/2) This model is 
equivalent to a 1-D Anderson impurity model. The inter- 
chain tunneling is associated with a transverse magnetic 
field acting on the impurity. For U > the magnetic sus- 
ceptibility and transverse magnetization (c^cr) will be 
rapidly enhanced with increasing U . Hence, the differen- 
tial tunneling conductance will be significantly increased 
by a large positive U. 

The above discussion of the exciton effect has been 
performed to leading order in |T| 2 . We will now go be- 
yond leading order in \T\ 2 and show that if |T| large and 
e > 0, that the Coulomb blockade is destroyed. To do 
this, we integrate out the particlc-holc excitations within 
the grains. This gives an effective action in imaginary 
time of the form 
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where a(r) = a ( ^^ff^pj j and a = h/(2ire 2 R T ). 

This is a one-dimensional XY model with long range 
interactions. This model had been studied within the 
framework of the renormalization group by Kosterlitz 
|l2f who found an order-disorder transition at ao — a c — 
2/(eir 2 ) for e > 0. The model is disordered for e < al- 
though the absence of an ordered phase when e = was 
a source of controversy. [|l3| — P~5|| . 

In order to understand the nature of the two phases, we 
calculated the conductance of the model using the Kubo 
formula G(ui) — (\I t (uj)\ 2 ) /u> where It(w) is the tunneling 
current. To leading order in an expansion in powers of 
1/a, a spin-wave calcualtion reveals that 



G = 



2V^a r((l + e)/2) / E Q 



Rq r(l + e/2) \nk B T 



(3) 



where E Q = e 2 /2C and R Q = h/e 2 = 4.11ML We see 
that G diverges at T — > 0. Hence one identifies the or- 
dered phase as subohmic. One can also calculate the 
conductance to leading order in ao. In this case one 

finds G ~ -^2- ( %^") 2(1 £) which vanishes as T -► 

indicating that the disordered phase is insulating. 

To understand the transition in greater detail, we have 
evaluated the DC conductance using a Monte Carlo sim- 
ulation fij). Using 111 



G = 



WRq 



h0 



7£ (t)(cos(0(t)-0(O))) dr (4) 



where 7 e (r) = [it^sT / Eq)/ sm(7rfcsrr?i)] _e , we evalu- 
ate G using (cos(</>(t) — <f>(r'))) obtained from a series of 
simulations including 10 5 cycle runs for the 10, 20, and 
32 timeslice systems and 2 x 10 5 cycle runs for the 64 and 
128 timeslice systems. The results for e = 0.2 are pre- 
sented in Fig. g The transition between the sub-Ohmic 
(high ao) phase to the insulating (low ao) phase is evi- 
dent. Interestingly one observed the conductance curves 
cross at a single point. This point identifies a transition 
at a critical value of a c = 0.9 which compares well to 
a c = 1.01 obtained by Kosterlitz RG treatment. 
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FIG. 3. Conductance obtained from a MC simulation. 
Notice, G(ao) curves intersect at G c ~ lle 2 /ft. This sug- 
gests that the conductance of the critical state is a universal, 
temperature independent constant. 

The fact that curves intersect at all indicates that 
a c separates metallic from insulating behavior. How- 
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ever, the observation that all line cross at G c w 11-^ « 
l/(2.3fcf2) seems to indicate that the critical state has a 
finite temperature independent conductance. The exis- 
tence of a critical state with a finite resistance can be 
understood as follows. In general, finite size scaling the- 
ory implies that the critical states exhibits correlations 
of the form (exp i<£(r) exp -#(0)) = (^.) d " 2+ '' F{hp/r) 
where tq — Ti/Eq is the width of the time slices, d = 1 
is the space-time dimensionality, and where F(x) is a 
universal scaling function which is finite as x — » 0. Ac- 
cording to Fisher, Ma and Nickel ]Ti| , r) = 1 + e is exact 
for our model. It follows that 



G c 



e 2 



(5) 



at criticality. Since the DC conductance depends only 
on the lo — > limit of the model, the above result is a 
universal but e dependent result. Typically simulations 
for other values of e reveal that 1/G C ~ 10 5 f2 except in 
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the e — » limit where G c — ► oo. It should be mentioned 
that the Herzog data reveals that the highest resistance 
metallic state in Au (width: 400A), Pbo.s5Bio.15 (width: 
575A and 850A), and Sn (width: 550A) wires have criti- 
cal resistances of 2/cO, 2kfl, 4fcf2, and 4fcO, respectively. 
Such values of G^ 1 are consistent with eqn. ||. 




FIG. 4. Modeling a crack as a series of parallel tunnel 
junctions (see inset). Observe that the charging energy de- 
creases with increasing wire width. The decreasing Eq causes 
the resistance gap to close and the resistance of non-critical 
subohmic wires to increase. 



Discussion - We wish to consider the relevance of our 
tunnel junction to the MI transition observed by Her- 
zog et al. 0| . Herzog's MI transition occurs in many 
different materials and exhibits a resistance gap which 
depends on the wire width, 500A — 6000A. This sug- 
gests that the transition does not involve some micro- 
scopic phenomenon in the interior/surface of the grains 
but, instead, involves some low energy collective behav- 
ior which requires two or more grains. In addition, it 
appears that the MI transition is not a percolation tran- 
sition in which the drop in the resistance is determined 
by purely geometrical considerations such as number of 
parallel connections spanning a backbone cluster. This is 
clear since percolation transition can occur in films but 
not in wires. At this point, one might model the wire as 
a disordered network of tunnel junctions where Coulomb 
blockade effects could produce an insulating phase. This 
is reasonable since charging energies are estimated Q 
to be ~ 100A". Assuming the tunnel junction conduc- 
tances to be broadly distributed, the wire resistance will 
be dominated by the tunnel junctions with highest resis- 
tance @. 

Now consider the dependence of the wire resistance on 



the wire width. Let 1$ be the phase coherence length i.e. 
the length scale that correlations (e^^e -4 ^ )) die off in 
the disordered phase. Then for wire widths w C i^, 
the low frequency behavior of the obstruction will be 
described by the tunnel junction model with a charg- 
ing energy which decreases with increasing wire width. 
For instance, in a narrow wire (w -C l<f>) phase differ- 
ence across multiple tunneling sites spanning a crack or 
weak-link will be equal. (See fig. 4.) Hence, the capaci- 
tances comprising the weak-link add in parallel. This is 
a useful observation since the decreasing charging energy 
associated with increasing wire width will cause the re- 
sistance of the non-critical subohmic states to increase 
like (T/EqY (eq. ||) similar to the behavior observed by 
Herzog et al.. The increase of the metallic wire resistance 
with wire width is a unique phenomenon which is difficult 
to obtain from alternative models. The model also pre- 
dicts that the resistance of insulating wires will decrease 
with increasing wire width. Consequently the resistance 
gap will decrease until w ~ 1$ where the crossover to two 
dimensional transport will occur. 

At this point, we should mention that our model makes 
several predictions which have not yet been tested. First, 
from the Kosterlitz RNG treatment demonstrates the 
existence of a characteristic energy scale, A = Eq{1 — 
oio/ac) 1 ^. This behavior should be detectable by exam- 
ining the differential conductance which obeys a scaling 
ansatz of the form dl/dV = ^(eV/E Q ) 5 F(eV/A) where 
5 = in order that scaling ansatz is consistent with eq. [|. 
The second prediction is that the transition to the phase 
coherent state in arrays of coupled tunnel junctions will 
be accompanied by a divergence in 1$. This unique be- 
havior should be readily observable in magnetoresistancc 
measurements performed on insulating granular wires. 

We should also mention that one should also be able to 
search for the subohmic to insulator transition in double 
quantum dot systems of the sort considered by Waugh 
et al. fll],^(J The double-dot systems have several useful 
features including (l.)a small number of tunneling chan- 
nels which implies the large phase shifts [0 required for 
an exciton effect, (2.) small intergrain capacitances |l9| , 
and (3.) precise control of the tunneling barrier between 
dots. Unfortunately, the Waugh experiment itself could 
not distinguish between a cross-over and a phase transi- 
tion. However, this is not an inherent limitation of the 
experimental method. So an attempt to search for this 
quantum phase transition in the double dot should be 
feasible and would certainly be most welcome. 
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